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I. 


INTRODUCTION 


In  1977,  Lloyd  and  Brown1  investigated  the  feasibility  of 
controlling  a  105mm  spinning  projectile  by  means  of  horizontal  and 
vertical  forces.  Their  numerical  calculations  yielded  the  surprising 
result  that  an  applied  constant-amplitude  yaw  moment  could  cause 
dynamic  instability.  This  result  was  then  explained  theoretically  by 
the  influence  of  a  nonlinear  term  in  <J>pp,  the  fixed-plane  coordinate 
system  spin  rate. 

The  Lloyd  and  Brown  work  was  limited  to  large  stability  factors, 
neglected  spin  and  density  variations  and  incompletely  considered 
gravity.  Reference  2  removed  these  limitations  and  thereby  gave  some¬ 
what  more  general  results  as  well  as  a  more  elegant  derivation.  Un¬ 
fortunately,  the  relations  of  Reference  2  were  applied  there  to  the 
case  of  a  nonspinning  statically  stable  missile.  Although  these  rela¬ 
tions  are  valid  for  statically  stable  missiles  with  spin  rates  that  are 
large  in  comparison  with  the  resonance  spin,  they  are  definitely  in¬ 
valid  for  zero  and  "near-zero"  spin.  It  is  the  purpose  of  this  memo¬ 
randum  report  to  give  a  more  rigorous  derivation  of  the  relations  of 
Reference  2  and  to  obtain  the  correct  equations  for  a  statically 
stable  missile  with  zero  spin. 


II.  GENERAL  THEORY 

Equation  (3.8)  of  Reference  2  can  be  written  in  the  form: 
l"  +  (H  -  i  P)  -  (M  +  i  PT)U  -  5e) 


=  ~  2  i  E[B//  +  (H  -  i  P)§']  (2.1) 


7.  K.H.  Lloyd  and  D.P.  Brown ,  "Influence  of  Gravity  and  Applied  Side 
Forces  on  the  Stability  of  a  Spinning  Projectile Weapons  Research 
Establishment  TR  1906 (W) }  South  Australia ,  November  1977 , 

AD  A053648.  (See  also  "Instability  of  Spinning  Projectiles  During 
Terminal  Guidance , "  Journal  of  Guidance  and  Control  2,  January- 
February  1979 ,  pp.  65-70.) 

2.  C.  11,  Murphy ,  "Effect  of  Horizontal  and  Vertical  Side  Forces  and 
Moments  on  Stability  of  a  Symmetric  Missile  in  Ascending  or  De¬ 
scending  Flight j"  USA  ARRADCOM  Ballistic  Research  Laboratory 
Memorandum  Report  No.  02925,  April  1979 ,  AD  A072308. 
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where 


E  =  (a/2Ke  tan0g 
£e  =  -  C/(M  +  i  P  T) 

a  =  [y  +  a  tan  B  ]  1 
6  e  e J 

and  the  other  symbols  are  defined  in  the  List  of  Symbols. 

The  right-hand  side  of  Equation  (2.1)  is  the  forcing  function  due 
to  <j>pp,  linearized  with  respect  to  the  equilibrium  angle  produced  by 

the  control  moment.  Since 

6  -  ee  =  d/2)  [i  -  le  *  l  -  le)  (2-2) 

§  can  be^eliminaced  from  Equation  (2.1).  The  result  has  a  quite  simple 

form  if  E  is  assumed  to  be  constant: 
e 

V  +  (H  -  i  P)E;  -  (1  -  i  E)(M  +  i  PT)(E  - 

=  -  i  Etc'7  +  (H  -  i  P)?'  ]  (2.3) 

where 

E  =  E/(l  +  i  E) 

% 

Equation  (2.3)  is  a  second-order  complex  equation  with  constant 
coefficients,  but  it.  involves  E  as  well  as  E  •  Equations  of  this  kind 
were  solved  in  Reference  3  by  assuming  a  solution  of  the  form: 

mi  s  ?  s  y  s  s 

E  =  Ee  +  r  +  k2e  2“  +  k4e  1  +  kge  2  (2.4) 


where  k.  and  4*.  are  complex  constants: 
1  3 


.3.  C,  H.  Murphy,  "Angular  Motion  of  Spinning  Almost-Syimetric  Missiles," 
Journal  of  Guidance  and  Control  2,  Ncve'mber-December  1979 ,  pp. 
504-510.  (See  also  USA  AEEADCOMballiatic  Research  laboratory 
Technical  Report  No.  02121 ,  November  1978 ,  AD  A0635Z8. ) 


If  Equation  (2.4)  is  substituted  in  Equation  (2.3)  and  the  co¬ 
efficients  of  the  four  exponentials  are  set  equal  to  zero,  there  results 


k.  ,  F.  +  i  M 

JULl  „  _J _ 2 _ 

kj  -  i[E  F.  ♦  Me] 


i  {  j  +  2  i  P(Vj  +  T)  ]  +  Me} 

F.  -  i  M  +  2  i  P(Y.  +  T) 
j  e  3 


(2.5) 


j  -  1,  2 


where  F.  =  +  (H-iP)'l'.  -  M-  iPT. 

3  3  3 


He  =  E(M  +  i  P  T) 


For  E  =  0,  =  kg  =  0  and  the  equation  F.  =  0  is  the  usual  quadratic 

equation  for  the  epicyclic  frequencies  and  damping  rates.  We  wish  to 
find  the  effect  of  small  nonzero  E.  Equation  (2.5)  can  be  written  in 
the  following  form,  where  4^  +  T  has  been  approximated  by  i<f>(  . 

ri(F.  -  2  P  *')  ♦  M  1  (E  F.  +  M  ) 

F.  +  i  Me  =  - 1 - i - ^ - 1 - —  (2.6) 

F.  -  2  P*'.  -  i  M 
3  3  e 

A  first  approximation  for  Fj  can  be  obtained  by  setting  the  small  right 
side  of  Equation  (2.6)  equal  to  zero. 


F.  =  -  i  M 
3  e 


A  better  approximation  for  F.  now  follows  from  Equation  (2.6)  by  re¬ 
placing  F.  on  the  right  side-’by  -  i  M  and  neglecting  F.?  teims  compared 
3  e 

to  E  terms. 
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For  this  approximation. 


b 


E  [2  P  <j>'  -  M  +  i  PT] 

i  + - - - — - 

2  P  +  i  (M  M  ) 
j  e  e 

-i  E  (2  P  -  M  -  i  PT) 

2  M'  -  i  ( H/  Me) 


(2.8) 


(2.9) 


III.  SOLUTION  FOR  SMALL  SPIN 


For  small  spin,  the  PT  term  in  Equations  (2.8  -  9)  can  be  omitted 
and  these  equations  reduce  to 


F.  =  E  M 


-l  + 


E  (1  -  Rj) 

1  +  i  (E  +  E)  Rj 


(3.1) 


k.  , 
J+3 

k. 

J 


-i  E  (1  -  R  ) 

1  -  i  (E  +  E)  Rj 


(3.2) 


where  =  M/(2  P  c|>(). 


The  iterative  process  that  produced  Equation  (3.1)  is  valid  only 
if  the  absolute  value  of  the  second  term  in  that  equation  is  smaller 
than  the  absolute  value  of  the  first  term.  For  small  |e|,  this  is  the 
case  if  ]R^|  is  less  than,  say,  unity  (a  conservative  upper  bound): 

I Rj I  <1  (3.3) 

This  condition  is  always  satisfied  for  the  rapidly  spinning  shell  of 
Reference  1,  but  its  applicability  to  a  slowly  spinning  f inner  must 
be  analyzed  by  use  of  the  concept  of  resonance  spin. 

Resonance  spin,  <(/  ,  is  defined  by  the  relations 

r  reson  ' 


re  son 


(3.4) 
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Relation  (3.3)  assumes  a  very  simple  form  in  terms  of  resonant  spin. 


/ 

reson 


> 


(3.5) 


Moment  of  inertia  ratios  can  be  as  high  as  10:1.  Thus,  relation  (3.5) 
would  be  satisfied  for  spins  in  excess  of  five  times  resonance. 


When  condition  (3.5)  is  satisfied,  the  magnitudes  of  k^  and  k?  are 

given  in  terms  of  k^  and  k2  by  Equation  (3.2)  and  are  small  for  small 

| £ | .  For  a  trim  of  12°  and  a  maximum  climb  or  descent  angle  of  45°, 

( E |  would  be  about  0.2.  For  these  conditions,  k4  and  k?  could  be 

neglected  and  the  motion  approximated  by  a  tricycle.  This  is  what  was 
done  in  Reference  2  and  Equation  (2.7)  yields  the  frequencies  and 
damping  rates  of  that  report. 

For  spin  rates  lower  than  five  times  resonant  spin,  Equation  (2.7) 
loses  validity  and  should  be  used  with  caution.  This  wau  not  done  in 
Reference  2  and  erroneous  results  were  given  for  zero  spin. 


IV.  SOLUTION  FOR  ZERO  SPIN 
For  P  =  0,  Equation  (2.5)  requires  that 

(Fj  -  b  M)Fj  =  0  (4.1) 

where 

b  =  (a  tan  6  )/v 
e  e'  1  e 

If  ag  is  not  zero.  Equation  (4.1)  has  two  distinct  roots: 

F.  =  0,  b  M  (4.2) 

The  four  roots  of  these  two  quadratic  equations  are  the  T1,  Y  Y  j,  Y  r 
of  Equation  (2.4).  If  we  identify  Y .  and  Y  ?  as  the  roots  of  the  first 
of  these  quadratic  equations  and  Y2  and  Y 2  as  roots  of  the  second,  then 

Fj  =  0,  F2  =  b  M  (4.3) 
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For  these  roots.  Equation  (2.5)  reduces  to 

k  =  -  k,  (4.4) 

4  1 

and 

k5  ■  -  '£A)S2  <4'5> 

Hence  the  solution  to  Equation  (2.3)  for  P  =  0  and  /  0  is 

i<i>  -1$,  T  i*,  .x  -i*, ”1 

C  -  +  K  (e  1  -  e  V  K  I  e  2  -  U Jt  )e  2|  (4.6) 

e  i  z  l  e  e  j 

where,  from  Equation  (4.3) 

Kj  “  Kjo  6Xp  C_Hs/2^ 
and  (ignoring  a  small  H2  term) 

*1  =  [-M]5* 

<  =  -  [-(1  +  b)M]^ 

For  the  special  case  a  =0  (which  occurs  when  the  fixed-plane 
transverse  control  moment  hts  no  horizontal  component) ,  Equation  (4.1) 
has  repeated  roots:  F  =  F2  =  0.  Equation  (2.4)  is  inadequate  to 

handle  this  situation.  Instead,  we  must  assume  a  solution  of  the  form 

¥  s  ¥  s 

5  =  Ee  +  (kL  +  k4s)e  1  +  (k2  +  kb  s)e  1  (4.7) 
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Direct  substitution  in  Equation  (2.3)  yields  the  relations: 


K  - 


k,  -  - 


/ms  ten  6 
e  e 


(k,  +  k2)  (4.8) 


Hence,  for  P  =  0  and  a  -  0,  vie  have 


i4>,  i#; 

C  =  3e  +  Kle  +  K2e 


1  M  8  tan  0 

e _ e_ 

2  v  <J>/ 

e  i 


(Kj  sim}>1  -K2sini{)2)  s 


(4.9) 


where 


>[  ■  -  *2  *  [  — M] 


We  see  that  in  both  zero-spin  solutions,  (4 . 6)  and  (4.9),  £  has 
no  effect  on  the  damping  rates  and  Pe  has  no  effect  on  the  frequencies. 
A  nonzero  &e  will  effect  just  one  of  the  frequencies;  zero  oe  will 
cause  secular  terms  to  appear  in  the  solution. 


V.  SUMMARY 

1.  The  relations  given  in  Reference  2  are  valid  for  gyroscopically 

stable  missiles  for  l<s  <<  °°  and  for  statically  stable  missiles  for 

8 


which 


reson 


>  5. 


l  i  tvji  jLQLicoii)  otaL'i'/  iiuooxiLC  nx  in  iu  wci  o^fx.n  i  a  xco ,  Ln^ 

motion  is  quite  complex  and  should  be  studied  by  numerical  simulation. 
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USER  EVALUATION  OF  REPORT 
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